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I.  INTRODUCTION 

This  report  is  the  second  annual  interim  report  on  contract  F44620-74-C-0057 
to  study  the  application  of  ultrasonic  acoustic  wave  scattering  to  the  detection, 
identification  and  characterization  of  defects  in  solids.  The  contract  was 
initiated  on  1 February  1974  has  entailed  approximately  one  man-year  of 
effort  equally  divided  between  the  theoretical  analysis  of  the  scattering 
process  and  an  experimental  program  to  verify  its  validity  and  measurability. 

II.  SUMMARY  OF  PROGRESS 

1.  Scattering  of  Acoustic  Waves  by  Spherical  Obstacles  and  Voids 

The  specific  goal  of  the  present  program  is  the  characterization  of  the 
scattering  of  acoustic  waves  by  ellipsoidal  obstacles.  However,  in  order 
to  verify  the  experimental  techniques  and  to  develop  the  necessary  theoreti- 
cal basis,  the  initial  task  chosen  was  the  review  of  the  existing  calculations 
of  scattering  by  spheres  and  the  verification  of  the  experimental  technique 
using  spherical  cavities  and  spherical  Inclusions  embedded  in  a solid  matrix.1 

The  results  of  this  work  have  been  partially  reported  In  a paper  sub- 
mitted to  the  "Journal  of  the  Acoustical  Society  of  America."  The  abstract 
of  this  paper  is  presented  here: 

a.  "Scattering  of  Longitudinal  Waves  Incident  on  a Spherical  Cavity  in  a 
Solid,"  B.  R.  Tlttmann  and  E.  Richard  Cohen 

Experimental  data  are  compared  with  theoretical  calculations  for  the 
scattering  of  ultrasonic  pulses  from  a single  spherical  cavity  embedded 
In  a metallic  solid.  Except  at  high  frequencies,  the  data  are  In 
reasonable  agreement  with  theory,  bearing  out  the  notable  features  pre- 
dicted for  the  angular  dependence  and  the  mode  conversion  between 
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longitudinal  and  transverse  waves.  At  high  frequencies  the  disagree- 
ment can  be  accounted  for  by  the  lack  of  monochromicity  in  the  ultra- 
sonic pulse.  The  experiments  were  typically  carried  out  in  the 
regime  where  the  wavelength  Is  approximately  equal  to  the  size  of  the 
spherical  cavity  (0.8  mn)  on  samples  of  diffusion-bonded  Ti-6%  A1-4JV. 
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2.  Scattering  by  a Spherical  Inclusion 

During  this  past  year  the  experimental  portion  of  the  program  has  con- 
centrated on  the  measurement  of  the  scattering  characteristics  of  the 
tungsten  carbide  inclusion.  Sample  preparation  was  carried  out  In  a manner 
similar  to  that  associated  with  the  preparation  of  the  800ym  spherical 
cavity  except  that  a commercially  obtained  tungsten  carbide  ball  with  a 
nominal  diameter  of  1/32"  (794ym)  was  inserted  into  the  cavity  at  the  center 
of  the  sample  prior  to  diffusion  bonding.  As  before,  the  sample  was  given 
a polygonal  geometry  with  14  flat  faces  each  just  under  0.5"  wide. 

The  experiments  for  incident  longitudinal  waves  were  carried  out  as 
before  with  a matched  pair  of  Panametrics  broadband  transducers.  The  three 
pairs  of  transducers  were  designated  as  2.25  MHz,  5 MHz  and  10  MHz,  but 
because  of  the  broadband  character,  with  a Q of  approximately  5,  waveform 
synthesis  was  required.2  The  2.25  MHz  transducers  showed  a response  charac- 
teristic with  a full  width  at  half  maximum  of  0.8  MHz.  The  peak  response 
of  the  5 MHz  transducers  was  at  4.5  MHz  with  a full  width  at  half  maximum 
of  1.9  MHz  and  for  the  10  MHz  transducers  the  peak  response  was  at  7.4  MHz 
with  a full  width  at  half  maximum  of  2.6  MHz  (Fig.  1).  On  the  other  hand, 
the  shear  wave  transducers  used  for  generating  and  receiving  transversely 
polarized  waves  were  high  Q devices  with  well-defined  monochromatic  response 
curves  which  were  within  1 to  1.5*  of  2.5  Wz,  5.0  MHz  and  10.0  MHz, 
respectively. 

The  results  for  incident  longitudinal  waves  are  shown  In  Figs.  2,  3,  and 
4.  The  angular  distribution  for  L-*t  scattering  was  calculated  for  2.25  MHz. 
The  synthesized  curve  using  the  measured  response  of  Fig.  1 was  essentially 
Identical  to  this  curve.  The  curve  marked  L-*T  was  calculated  for  a frequency 
of  2.5  MHz  since  this  Is  the  frequency  to  which  the  receiver  responds.  In 
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Fig.  3 the  mode-converted  angular  distribution  Is  similarly  calculated  for 
a frequency  of  5.0  MHz,  while  the  direct  scattered  curve  (L-*L)  is  calculated 
from  the  spectral  synthesis  (solid  curve).  Also  shown  are  monochromatic 
4 and  5 MHz  curves.  These  indicate  the  sensitivity  of  the  angular  distri- 
bution to  frequency  or  to  size  of  the  scatterer  as  well  as  the  averaging 
effect  of  frequency  synthesis.  In  Fig.  4 the  10  MHz  monochromatic  curve 
shows  a great  deal  of  structure,  in  particular,  the  surprisingly  sharp  dip 
at  92°  in  the  mode-converted  scattered  shear  wave.  The  directly-scattered 
10  MHz  curve  (M.)  is  indicated  although  it  has  little  bearing  to  the 
experimental  data  since  there  is  little  10  MHz  content  in  the  pulse.  The 
7.4  MHz  monochromatic  curve  (corresponding  to  the  peak  frequency  of  the 
broadband  pulse)  is  much  closer  to  the  synthesized  curve  which  shows,  as  is 
to  be  expected,  rather  broad  features  with  no  sharp  detail. 

The  experimental  data  was  fitted  to  these  curves  by  choosing  an  arbitrary 
normalization  in  order  to  give  the  best  fit.  The  lack  of  agreement  between 
the  synthesized  calculation  and  the  data  for  L-*L  scattering  in  Fig.  4 may 
be  the  result  of  scatter  from  a secondary  metallurgical  phase  which  would 
become  more  important  at  higher  frequencies.  One  possible  source  of  error 
in  the  mode- con  verted  cases  (L-*T)  is  the  unfortunate  mismatch  between  the 
transmitter  frequency  which  peaks  at  7.4  MHz  and  has  a distribution  with  a 
mean  frequency  of  7.26  MHz  with  a standard  deviation  of  1.16  MHz;  there  is 
thus  very  little  power  at  10  MHz  for  the  receiver  to  respond  to. 

Experiments  with  incident  shear  waves  were  carried  out  with  shear-cut 
quartz  plates  which  were  bonded  to  the  faces  of  the  sample  with  low-melting- 
point  wax  and  which  could  be  polarized  either  parallel  or  perpendicular  to 
the  scattering  plane.  The  results  are  presented  in  Figs.  5,  6 and  7.  The 
angular  distribution  for  the  mode-converted  scattered  compresslonal  wave 
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Fig.  3 the  mode-converted  angular  distribution  Is  similarly  calculated  for 
a frequency  of  5.0  MHz,  while  the  direct  scattered  curve  (L-*l)  is  calculated 
from  the  spectral  synthesis  (solid  curve).  Also  shown  are  monochromatic 
4 and  5 MHz  curves.  These  indicate  the  sensitivity  of  the  angular  distri- 
bution to  frequency  or  to  size  of  the  scatterer  as  well  as  the  averaging 
effect  of  frequency  synthesis.  In  Fig.  4 the  10  MHz  monochromatic  curve 
shows  a great  deal  of  structure,  in  particular,  the  surprisingly  sharp  dip 
at  92°  In  the  mode-converted  scattered  shear  wave.  The  directly-scattered 
10  MHz  curve  (L-M.)  is  indicated  although  it  has  little  bearing  to  the 
experimental  data  since  there  is  little  10  MHz  content  in  the  pulse.  The 
7.4  MHz  monochromatic  curve  (corresponding  to  the  peak  frequency  of  the 
broadband  pulse)  is  much  closer  to  the  synthesized  curve  which  shows,  as  Is 
to  be  expected,  rather  broad  features  with  no  sharp  detail. 

The  experimental  data  was  fitted  to  these  curves  by  choosing  an  arbitrary 
normalization  In  order  to  give  the  best  fit.  The  lack  of  agreement  between 
the  synthesized  calculation  and  the  data  for  L-*L  scattering  in  Fig.  4 may 
be  the  result  of  scatter  from  a secondary  metallurgical  phase  which  would 
become  more  important  at  higher  frequencies.  One  possible  source  of  error 
In  the  mode-converted  cases  (L-*-T)  is  the  unfortunate  mismatch  between  the 
transmitter  frequency  which  peaks  at  7.4  MHz  and  has  a distribution  with  a 
mean  frequency  of  7.26  MHz  with  a standard  deviation  of  1.16  MHz;  there  is 
thus  very  little  power  at  10  MHz  for  the  receiver  to  respond  to. 

Experiments  with  Incident  shear  waves  were  carried  out  with  shear-cut 
quartz  plates  which  were  bonded  to  the  faces  of  the  sample  with  low-melting- 
polnt  wax  and  which  could  be  polarized  either  parallel  or  perpendicular  to 
the  scattering  plane.  The  results  are  presented  In  Figs.  5,  6 and  7.  The 
angular  distribution  for  the  mode-converted  scattered  congressional  wave 
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(T-H.)  is  the  same  as  the  mode-converted  scattered  shear  wave  from  incident 
congressional  waves  (L+T).  This  "reciprocity"  in  the  scattering  is 
discussed  more  fully  in  Section  II. 3 below. 

The  curves  marked  T-*T|j  and  T-*T^  in  Figs.  5,  6 and  7 represent  the 
amplitudes  of  the  scattered  shear  wave  when  the  plane  of  polarization  of 
the  incident  shear  wave  Is  parallel  and  perpendicular  to  the  scattering 
plane.  The  measurements  were  made  in  these  two  orientations.  In  general, 
when  the  scattering  plane  Is  oriented  at  an  angle  <t>  to  the  plane  of 
incident  polarization,  the  scattered  shear  wave  is  elliptically  polarized; 
the  amplitude  of  the  shear  wave  in  the  scattering  plane  is 

A|j  (R,8,$)  = ^ (0)cos<J>  (2.1) 

and  the  amplitude  of  the  shear  wave  perpendicular  to  that  plane  is 

Ai(R,0,<j>)  3 | Ay^.  (6)s1n+  . (2.2) 

The  measurements  at  5 MHz  however  showed  the  cross-polarized  wave  at  per- 
pendicular scattering  to  have  an  appreciable  amplitude.  This  Is  not  under- 
stood at  this  time;  it  may  be  the  result  of  scattering  by  the  metallurgical 
second  phase  In  the  titanium  alloy.  This  background  scattering  may  also  be 
the  cause  of  the  lack  of  qualitative  agreement  between  theory  and  experiment 
observed  In  the  shear  wave  scattering  at  5 MHz,  and  requires  further 
investigation. 
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3.  Reciprocity 

In  the  course  of  the  measurements  of  mode-conversion  in  scattering 
from  spherical  cavities  and  inclusions,  it  was  recognized  that  the  angular 
distribution  for  the  mode-converted  component  is  the  same  for  the  shear 
wave  scattered  from  an  incident  compressional  wave  and  for  the  compressional 
wave  scattered  from  an  incident  shear  wave.  Thus,  in  Figs.  5,  6 and  7,  the 
curves  labeled  T-*t  differ  only  by  a constant  factor  from  those  in  2,  3 and 
4 labeled  L-*-T. 

This  result  had  not  been  anticipated  but  is  contained  in  the  analytic 
expressions  for  the  scattering  amplitudes.  The  amplitudes  of  the  scattering 
by  a spherical  scatterer  are  given  by  expressions  of  the  following  form: 

For  incident  compressional  wave,  normalized  to  unit  amplitude  for  the 
displacement  in  the  incident  wave3: 

ur  ■ -ef-  -y^(2n-H)AHPrl(cose)  (3.1) 

^^2J(2"+,)BnVcosS)  (3-2) 


For  incident  shear  waves,  similarly  normalized: 


ie 


-ikr 


ir^E&Srr  v;<—> 


1e'krfn_(fcV2n+1 

C0StL»n TH+TT 


ie 


<r 
-iter 


p!(cos9)  dP1 (cos8) 
u n . p n 

*n  sine  rn  de 


4 icr 


‘,n£s7snT 


.1 
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dP  (cose)  P‘(cos9) 

( D + r Jl3 

fl  d8  n sine 


(3.3) 

(3.4) 

(3.5) 
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By  direct  evaluation  of  the  equations  which  define  the  constants 
in  these  expansions3  we  can  show  that 

En  - n(n+l)(k/<)3Bn  . (3.6) 

This  result  can  be  obtained  without  explicitly  evaluating  the  complete 
expressions  for  En  and  Bn-  It  had  previously  been  shown  to  hold  for  the 
Rayleigh  limit,  ka  « 1,  but  is  now  established  for  arbitrary  values  of  ka. 
Thus  we  have 

u(shear)/uQ(compr. ) = u(T-»l)/u(L-*-T)  = k2/<2  (3.7) 

r y 

and  the  differential  scattering  cross  sections  0^(9)  and  0^.(9)  are  also 
in  the  same  ratio 

aT-L(9>'VT<6>  * k^2  ' ' (3'8) 

This  result  in  fact  appears  to  be  a special  case  of  a much  more  general 

reciprocity  theorem  valid  for  a scatterer  of  arbitrary  shape  in  an 

elastic  medium,  which  follows  from  general  arguments  of  symmetry  and 

invariance.  If  we  consider  an  arbitrary  shape  (Fig.  8)  with  an  incident 

plane  wave  of  polarization  p and  propagation  vector  Kj,  the  scattered  wave 

amplitude  with  polarization  q and  propagation  vector  can  be  written  as 
2 

Kjfpq(KpKj)  and  reciprocity  is  contained  in  the  invariance  of  f to  the 
interchange  of  p and  simultaneously  Kj  <+  -K^,  i.e., 

VW  ■ V'V-V  <3-9> 
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For  the  sphere  the  additional  symmetry  is  such  that  f ( K.j ,K^ ) depends  only 
on  the  angle  of  scattering;  i.e.,  on  K.^ - K.^  and  not  on  the  orientation  of 
the  vectors  individually. 

The  reciprocity  theorem  in  its  most  general  form  may  easily  be 

established  for  an  arbitrary  elastic  medium. u We  consider  a region  of  space 

with  density  p and  elastic  coefficients  C (both  functions  in  general 

of  position).  The  elastic  displacements  at  point  r.  produced  by  a unit 

delta-function  fc  ce  of  frequency  co  applied  in  the  j-direction  at  point  r^ 

★ 

is  a solution  of  the  equation 


(Cmnps(r)us,p(^l>)  .n  * * V*--!1  ’ ° (3',0) 

Similarly  we  write 

(W(r)us,p(M<!>)  ,n  * = ° <3-'” 


Now  multiply  the  first  equation  by  uq(r;r„),  the  second  by  -uJ(r;r, ),  add 

m c Ri  i 

and  integrate  over  the  region  (R.  Because  of  the  symmetry  of  the  elastic 

coefficients,  C(r)  = Cn.Ar)  = C m„„(r),  one  obtains 
mnps  — psmn ' nmps  — 


mnps 


-"i(D!:,)Clmlps^,p(r;r2)]dSn 


q 

S,P 


■ * u^r,> 


(3.12) 


★ 

We  utilize  the  ordinary  tensor  notation  In  which  summation  is  assumed  on 
repeated  indices  and  a comma  Indicates  differentiation  with  respect  to  the 
following  coordinate  indices. 
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If  the  region  fl  extends  to  infinity  the  displacements  fall  to  zero  at 
least  as  fast  as  1/r  and  the  surface  integral  vanishes.  If  the  region  fl 
is  finite  we  assume  homogeneous  boundary  conditions  such  that  either  the 
displacements  or  the  normal  derivatives  of  the  stresses  vanish  on  the 
boundary  surface  and  again  the  surface  integral  will  vanish.  We  therefore 
obtain  the  reciprocity  rule: 

u j (ht  ‘»r.2)  = (3.13) 

which  states  that  the  j -displacement  at  r^  due  to  a unit  q-force  at  r^  is 
equal  to  the  q-displacement  at  r^  resulting  from  a j-force  at  r 1 . Thus  a 
receiver  at  r^  responding  to  waves  with  polarization  q from  a transmitter 
at  r..j  exciting  displacements  of  polarization  p will  record  the  same  signal 
when  the  transmitter  and  receiver  are  interchanged. 
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4.  Integral  Formulation  of  Scattering 

As  an  alternative  to  the  normal  mode  expansion  which  is  so  useful  in 
spherical  coordinates  and  whose  extension  to  spheroidal  coordinates  involves 
serious  computational  problems,  an  integral  formulation  has  also  been 
explored,  based  on  the  utilization  of  Green's  function  for  an  elastic  medium. 

We  consider  an  infinite  medium  of  density  p'(r ) and  elastic  constants 
Cmnps(-)  such  for  r “*  00  the  medium  is  homogeneous,  isotropic  and  uni- 
form with  density  p and  elastic  constants  « X«nf1«ps  + «msSnp)- 

The  elastic  displacements  with  frequency  <d  in  the  medium  are  given  by 

f4-1) 

We  introduce  the  Green's  function  G^(r;r^)  for  the  uniform  medium  which 
satisfies 


Cmnp$®s,pr/— 1 ^ 4 ' 0 


The  inhomogeneity  of  the  medium  is  then  defined  by 


(4.2) 


ACmnps^ 


^mnps^  " Cmnps 


Ap(r)  * p'(r)  - p 


(4.3) 

(4.4) 


In  a manipulation  similar  to  that  in  the  previous  section  we  obtain 


SC579.3IR 


fa  }Gi(l’:il)(Cnnps(j:>us,P(-)),n  - umWCmnpsGs,pp(^l 


+ Ap(r)G^l(r;r1)um(r)jdv 


I 0 r,£fl 


(4.5) 


Then,  for  r^efl  one  obtains 


uj(rl>  ’fs  [c'mps(i)Gi(^^)us.p(t>  ' Wu»(-)Gs.P(-;!;l)]dSn 
+fj^p(r)ui2G^(r;r1  )um(r) 


“ACmnps(-)<^s,p(-;-l  )um,n^pv 


(4.6) 


We  wish  to  define  a situation  with  an  incident  unperturbed  wave  traveling 
through  the  uniform  isotropic  medium  characterized  by  the  constant  proper- 
ties p and  The  Inhomogeneity  (the  scattering  center)  characterized 

by  Ap  and  AC„^„„  are  limited  to  be  contained  in  a volume  Vot  the  volume  of 
mnps  Z 

the  scatterer.  The  Green's  function  for  the  uniform  medium,  given  by  Eq. 
(4.2)  may  be  shown  to  vanish  as  lr.-r.jl~1  as  r 1 -*00. . On  the  infinite  surface 
S,  far  from  the  scattering  center,  the  displacement  um(r)  is  just  the  inci- 
dent wave  ui°)(r).  Then 
m — 

uJ0)<n,>  * W { sJ (n-.r, >U<»> (r)  - u<°>(r)GJip(r;r,)jdSn  (4.7) 


and  hence 
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uj (—1 ) ’ *L  laoCr^G^tr-.r^u^fr) 

V2  1 

-‘‘mnpsWi.p^K.nWfc  ' (4'8) 

In  an  infinite  uniform  medium  G^(r;jr^)  can  be  a function  only  of  the 

difference  jr-r_-j  (and  if  the  medium  is  also  isotropic,  only  of  the  scalar 
distance  R = Ir-j^l)-  Thus,  we  can  rewrite  Eq.  (3.13)  as 

eiter,)-«;(£,ir)-  VR)*  VR)  (4'9) 


where  R = | • Equation  (4.8)  can  then  be  written 

uj(r)  * Uj°^(r)  + y j u2Ap(rOGjm(r-r^um(r') 


+ ACmnps(^Gjm,n^up,s(^[d^ 


It  may  shown  from  Eq.  (4.2),  for  the  Isotropic  medium5 


VR) 


[*  SJ"  * 3xJ3*m  ) 


.2  . -iieR  .2  -IkR 

3 \ e 3 e 


R 3x.3x„  R 
j m 


with  pw2  » k2u  * k2(X+2u). 


(4.11) 


From  the  solution  of  Eq.  (4.8)  one  may  determine  the  scattered 
amplitudes  by  writing 


Uj(r)  ■ Uj°^(r)  + ujs^(r) 


SC579.3IR 


and  the  asymptotic  expression  for  Gjm(R)  for  kR  » 1.  The  scattered  wave 
amplitude,  written  for  kR  >>  1 , (i.e.,  the  far-field  limit)  can  be  expressed 
as 


(4.13) 


in  which  U^(fi)  is  the  amplitude  of  the  longitudinal  wave  and  U^(ft)  that 

of  the  transverse  wave.  These  amplitudes  are  functions  of  the  scattering 

direction  fl,  .(R  = fiR),  as  well  as  the  direction  and  polarization  of  the 

incident  wave  u^(r). 

m — 


(r)Ap(r)e1*<— d£ 


♦•[»  yv  j 

Tx^bk  f Cuj,j(!l)4X(l)  + 2iHnj£»p«j,p(n)] 


(4.14) 


and 


(4.15) 
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5.  The  Born  Approximation  and  the  Rayleigh  Approximation 

One  can  write  equation  (4.10)  and  its  derivative  in  a form  which  gives 
a set  of  coupled  integral  equations  for  the  displacements  and  the  elastic 
strains,  a ^{u^  j + Uj  j).  ^ have 

uj(l)  a ^o)(r)  + w2y*Ap(rOGjm(r-r^um(rOdr' 

+ /aC™ npSSj!l<n-r3eps(r-)<Jr-  (5.1) 

where  Gjl^n(R.)  is  the  symmetrized  Green's  function  derivative. 


Gj!i,(R>  * ? 


G,  (R)  + G.  (R) 
jm,n  — ' jn,m'— ' 


Equation  (5.1)  is  equivalent  to  (4.10)  because  of  the  symmetry  of  Cmnp$. 

From  this  equation  one  can  now  write  an  equation  for  by  differentiation 

but  we  must  Include  the  contribution  of  the  delta-function  In  gP^  . To 

j ;mn  ,s 

simplify  the  notation  we  now  introduce  the  six-dimensional  vector  representa- 
tion for  the  strains.  (The  components  of  the  six-dimensional  vectors  will  be 
designated  with  Greek  subscripts.)  We  also  Introduce,  for  convenience,  the 
dilatation 


D(I>  ■ Vm<H>  5 el<t)  * *2(r)  * s3(r) 

and  write 


(5.2) 
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Uj(r)  * uj°^(r)  + i/J'tip ( l) G ( r- 0 um ( r) dr 
+ yAX(rOGj°^r-r')0(r')dr' 

+ 2 ^(rOG^Cr-r^egCOdr'  (5.3a) 


Differentiating  this  equation  and  symmetrizing  the  resultant  expression 
gives  us 


:a(r)  -e^(r)  + oj^Ap ( rf)  G^ ( r-rO um ( r1) dr' 
+J&\  ( r1)  G^  ^ ( r-  r)  D ( r)  dr' 

+ zJlm  (r!)G^  ^ (n-r-Vg  (Odr*  . 


(5.3b) 


where 


:(o) 


g(2) 

can 


(R>  * WR) 

(5.4a) 

<R>  * si!>>  • 

r 

(5.4b) 

<R>  ■ Sjp;mtR)  ’ l|S,,p<R>  * WR) 

■ G";jp(R> 

(5.4c) 

(R)  - }[Gj(°>iR)  * G<°j(R)] 

(5.4d) 

(R)  • S<J>8(R)  ■ }[«]’»„«>  ♦ $>j(R) 

(5.4e) 

If  we  use  the  Incident  wave  field  uj0^  In  equations  (4.14)  and  (4.15) 
we  obtain  the  Born  approximation  as  It  has  been  developed  by  Krumhansl , 
et  al.' 
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[a]  V* 

For  incident  longitudinal  wave  uj  ' * -j*- e where  K is  the  incident 

2 

wave  vector,  KjKj  * k . Then  the  scattering  angle,  9 is  given  by 
fijKj  * kcose  and  we  have 


,U)  „ k_  A£ 


♦ 2Aucos2el  /j(kn-K)-r 
x + 2y  j ye  aL 


(5.5a) 


}jen<a-K)-r 


(5.5b) 


For  an  incident  transverse  wave,  Uj0^  3 u^e"*— — with  | 
and  K'jiij  * 0,  and  * sin8cos$.  Then 

u<*>  ■ £ sinecos*  . 2 t & cose]/^-^  dr 


r.  UjUj  • 1 . 


(5.6a) 


4°  —K  * ^ dr 


(5.6b) 


The  components  parallel  and  perpendicular  to  the  scattering  plane  may 
be  found  by  projection  of  the  vector  U^,  and  one  obtains 


-■f-y - cos$  £2-  COS0  * ^ cos29  /e1  dr 

4ttR  [ p y J — 


(5.7a) 


u^s^(t*tA)  » - 


2 -1<R 


sin*  tot&cose]/e1(^'idr 


(5.7b) 


This  formulation  of  the  Born  approximation  Is  adequate  only  for  the  case 
where  |dp|  « p,  |AX|  « X and  |dy|  « y,  otherwise  the  delta  function 
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contributions  from  and  G^  in  equation  (5.3b)  cannot  be  ignored. 

In  order  to  identify  these  delta  functions  we  need  look  at  the  Green's 
function  in  the  limit  <R  « 1 and  we  write 


V*> 


(5.8) 


where  the  omitted  terms  are  finite  as  R goes  to  zero.  To  the  same  order  of 
approximation  we  can  now  evaluate  the  second  order  derivatives,  but  to 

distinguish  the  delta  function  we  introduce  convergence  parameter  a 'nd  let 

2 2222222 
R be  given  by  R * + a * r + a . After  some  straightforward 

algebraic  manipulation  one  finds 


Gjm,ps^ 


a 

T7 

8irpo)  R . 


a2(<2  - k2)&.  - 2k2R25.  5 + 

jmps  jm  ps 


(5.9) 


where  Ajmps  is  the  comPletely  symmetric  tensor  djmps  - 6jm5p$  + 5jp6m$  * 

6-5 and  in  which  we  have  also  omitted  terms  which  are  spherically  ani- 

js  pm 

sotropic  around  the  point  R a 0 and  for  which  the  Integration  over  angles 
vanishes.  Except  for  the  point  r * o it  is  obvious  that  equation  (5.9) 
vanishes  as  a goes  to  zero.  The  delta  function  is  then  extracted  by 
writing 


Gjm,ps^ 


«(B.)  r 

2pu 


W‘2  • k2)Jl  - 2‘2VpSJ2  + 


(5.10) 


where 
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I 


f ' r*~dr  _ r°  u2du  _ _2_ 

Jo  IT -Jo  (1  ♦ U2)7/2  ■ 15 


J,  • 11m  a4  fr 

a -m3  Jo  R 


, _ ,a_  -2  r r2dr  _ /*"  u2du  _ 1 

**  "'JIT  Jo  1777?^  3 


We  then  obtain 


S(R) 


(<2  - k2)A^mnc  - 5^{,Jnc 
jmps  jm  ps 


and 


(5.11) 

(5.11) 


GjP»Pm 


6jm6(R)  + . . . 


(5.11a) 


G . 

jm,  pp 


G j p > pm 


2k4 


3po> 


V®  * • • 


(5.11b) 


These  expressions  satisfy  the  delta-function  component  «f  equation 
(4.2).  We  note  that  the  omitted  terms  in  equation  (5.11)  may  and  indeed  do 
contain  singularities  of  order  R"^  but  these  lead  to  finite  contributions 
to  the  integrals  in  (5.3b).  The  si ngif leant  feature  of  the  delta-function 
contributions  in  equation  (5.3b)  is  that  they  are  independent  of  the 
geometry  of  the  scatterer.  Equation  (5.3b)  may  then  be  written  in  the  form 


ea(l)  » e^(r)  - Aageg(r)  + J J"  Ap(rO  G^(  |r-rj|  Ju^rOdr' 

+ -ft A(r}fi^3) ( j r-rT  )D(r)dr' 

♦ 2 yAu(»36^)(|r-r3)e6(r5dW  (5.12) 
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-(2) 

where  A _(r)  is  the  contribution  of  the  delta  functions  and  Gv  '(R)  and 
eng  — a 

G^(R)  are  the  regular  (i.e.,  delta-function-free)  parts  of  the  Green's 

ctp 

function.  We  now  let  S denote  the  "shielding"  matrix  which  is  the  inverse 

of  5 , + A fl(r)  and  write 
a8  aS  — 

ea(r.)  ■ Say(r)  e^(r)  + J" ( r-r  )u|n(r')dr' 

+ )D(r')dr.> 

+ 2y"Ay(r)g^)(|r-r3)eB(r)dr;  (5.12a) 

Our  "corrected"  Born  approximation  involves  dropping  the  contributions  of 
the  integrals  in  this  equation  and  using  as  our  first  approximation  the 
"shielded"  strain  field 

e (r)  * S (r)e(o)(r)  (5.12b) 

a — ay  Y “ # 

in  equations  (4.14a-b)  to  calculate  the  scattered  amplitudes. 

There  is  in  fact  no  need  to  evaluate  the  matrix  S(r);  it  has  already 
been  evaluated  in  the  exact  solution  of  the  spherical  scatterer.  In 
lowest  order  the  shielding  matrix  can  be  extracted  from  the  Rayleigh  limit 
in  the  spherical  problem.  Since  S(r)  is  geometry-independent,  the  Rayleigh 
limit  can  then  be  converted  Into  the  Born  approximation  when  one  expresses 
the  scattering  amplitudes  in  the  form 
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It  is  clear  that  a3/ 3 for  the  sphere  is  simply  V/4-rr  and  hence  should  be 
replaced  by  the  appropriate  Fourier  transform  integral;  e.g., 

^/3  -*  ^ - dr  in  equation  (5.5a). 

We  have  omitted  some  terms  in  evaluating  the  delta-function  contributions 
in  equations  (5.3a,b)  arising  from  the  spatial  derivatives  of  the  elastic 

constants,  which,  even  for  a scattering  volume  of  constant  properties,  will 
not  vanish  because  of  the  step  jump  in  properties  at  the  surface.  This 
will  introduce  a surface  integral  into  the  shielding  matrix  A^g  in  equation 
(5.12).  Such  a term  is  to  be  expected  because  the  shielding  must  depend 
on  both  the  shape  and  the  orientation  of  the  scattering  object.  Such 
effects  are  well  known  in  both  the  static  elasticity  problem  and  in  the 
analogous  scattering  of  electromagnetic  radiation.  For  a spherical  scatterer 
these  terms  can  be  exactly  included  by  using  the  Rayleigh  limit  of  the 
exact  calculation  to  evaluate  the  shielding  factor.  For  non-spherical 
scatterers,  however,  equation  (5.12a)  is  incomplete  and  is  appropriate  only 
for  shapes  which  differ  only  slightly  from  spherical.  It  is  expected  that 
the  complete  identification  of  all  of  the  delta-functions  in  equations 
(5.3a)  and  (5.3b)  will  yield  a significantly  improved  formulation  that  will 
allow  the  Born  approximation  to  provide  an  appropriate  first  approximation 
to  scattering  by  such  extremely  non-spherical  defects  as  cracks  or  discs. 

Using  the  Rayleigh  limit  for  the  exact  solution  to  the  problem  of  the 
spherical  scatterer  allows  us  to  define  a surprisingly  simple  correction 
to  the  8orn  approximation. 
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We  define  the  coefficients 


A = 


AX  + | Au 


0 X + 2u  + AX  + y Ay 


A2  = 


Ai 

w + (3  + 2k  /<  )Ap 


We  then  obtain,  for  the  scattered  wave  amplitudes  with  incident 
compressional  waves 


u^(iWO 


kVikR 

4ttR 


A£ 


cos9  - A.  - 


— 2~  A2(cos^0 
< 


uyei(ka-K)-r  ^ 

(5.13a) 


/ \ 2 “ i kR 


^ ^ A2cose|sin9 yV^-y-tdr  (5.13b) 


We  note  that  in  this  approximation  the  scattered  transverse  wave  is 
polarized  in  the  scattering  plane.  For  a scatterer  with  axial  symmetry 
around  the  axis  of  the  incident  wave  vector  this  would  follow  from  simple 
symmetry  arguments.  The  implication  here  is  that  for  an  arbitrary  shape 
the  out-of-plane  polarization  amplitude  is  geometry-depencent  and  hence  of 
order  < V e (where  e is  a measure  of  the  asymmetry)  smaller  than  the 
in-plane  amplitude. 

For  the  scattered  wave  amplitudes  with  incident  shear  waves  we  obtain 


t2a-1kR 

u'  '(t-*-£)  * 4^ — sinecostf 


.(*) 


2 -i<R 
< e 


A£ 

p 


*k-'-*'-dr  (5.14a) 


u (t~V  ■— 4^-cos» 


2 -IkR 

<*  (t*ti>  ■ * far  s1n* 


- f a2 cosejyi 1 

cose  - A2(2cos20  - dr 


^-A2cose 


(5.14b) 


(5.14c) 
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The  comparison  of  the  first  order  Born  approximation  (equations  (5.5a 
(5.6a),  and  (5 . 7a , b ) ) and  the  corrected  Born  approximation  (equations 
(5.13a,b)  and  (5.14a,b,c))  with  the  exact  calculations  for  the  scattering 
from  a WC  inclusion  in  Ti  alloy  is  shown  in  Figs.  9-12.  It  is  clear 
that  in  this  case  in  which  we  have  Ap/p  = 2.1,  AA/X  = 0.8  and  Ay/y  = 4.4 
the  "shielding"  effect  should  be  quite  important.  As  expected  for  a 
situation  with  <a  not  small,  (<a  = 1.87)  the  Born  approximation  is  not  too 
good,  particularly  in  the  forward  scattering  direction.  The  corrected 
Born  approximation  is  significantly  better. 
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Fig.  1 

Fig.  2 
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Figure  Captions 

Frequency  analysis  of  the  transmitted  pulses  using  nominal  2.25  MHz, 
5.0  MHz  and  10  MHz  broadband  transducers. 

Scattering  distributions  for  incident  compressional  waves;  2.25  MHz. 
Scattering  distributions  for  incident  compressional  waves;  5.0  MHz. 

Scattering  distributions  for  incident  compressional  waves;  10  MHz. 

Scattering  distributions  for  incident  shear  waves;  2.5  MHz. 

Scattering  distributions  for  incident  shear  waves;  5.0  MHz. 

Scattering  distributions  for  incident  shear  waves;  10.0  MHz. 

Reciprocity  in  scattering  of  acoustic  waves. 

Comparison  of  Born  approximations:  Direct  scattering.  Incident  com- 
pressional wave,  ka  = 0.89,  WC  sphere  in  Ti64.  , exact  solution; 

, lowest  order  Born  approximation;  — shielded  Born  approxima- 
tion. 

Comparison  of  Born  approximations:  Mode  conversion,  longitudinal  to 

transverse  or  vice  versa,  ka  = 0.89;  <a  3 1.87;  WC  sphere  in  T i 64 . 

, exact  solution;  lowest  order  Born  approximation; 

— , shielded  Born  approximation. 

Comparison  of  Born  approximations:  Parallel  scattered  shear  wave, 

<a  « 1.87;  WC  sphere  In  T164.  , exact  solution;  , lowest  order 

Born  approximation;  — , shielded  Born  approximation. 

Comparison  of  Bom  approximations:  Perpendicularly  scattered  shear 

wave,  tea  * 1.87;  WC  sphere  in  Ti64.  , exact  solution;  , lowest 

order  Born  approximation;  ---,  shielded  Born  approximation. 
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Fig.  2 
Fig.  3 
Fig.  4 
Fig.  5 
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Fig.  10 


Fig.  11 
Fig.  12 


Figure  Captions 

Frequency  analysis  of  the  transmitted  pulses  using  nominal  2.25  MHz, 
5.0  MHz  and  10  MHz  broadband  transducers. 

Scattering  distributions  for  Incident  compress ional  waves;  2.25  MHz. 
Scattering  distributions  for  incident  compress ional  waves;  5.0  MHz. 
Scattering  distributions  for  incident  compressional  waves;  10  MHz. 
Scattering  distributions  for  incident  shear  waves;  2.5  MHz. 

Scattering  distributions  for  Incident  shear  waves;  5.0  MHz. 

Scattering  distributions  for  incident  shear  waves;  10.0  MHz. 
Reciprocity  in  scattering  of  acoustic  waves. 

Comparison  of  Born  approximations:  Direct  scattering.  Incident  com- 
pressional wave,  ka  = 0.89,  WC  sphere  in  Ti64.  , exact  solution; 

, lowest  order  Born  approximation;  — , shielded  Born  approxima- 
tion. 

Comparison  of  Born  approximations:  Mode  conversion,  longitudinal  to 

transverse  or  vice  versa,  ka  * 0.89;  <a  * 1.87;  WC  sphere  in  Ti64. 

, exact  solution;  lowest  order  Born  approximation; 

— , shielded  Born  approximation. 

Comparison  of  Born  approximations:  Parallel  scattered  shear  wave, 

ica  « 1.87;  WC  sphere  in  T164.  , exact  solution;  , lowest  order 

Born  approximation;  — , shielded  Born  approximation. 

Comparison  of  Bom  approximations:  Perpendicularly  scattered  shear 

wave,  xa  * 1.87;  WC  sphere  in  T164.  , exact  solution;  , lowest 

order  Bom  approximation;  — , shielded  Born  approximation. 
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